Superconductivity and magnetic order in CeRhlns; spectra of coexistence. 



O 

o 

(N 

u . 

Ph. 
< 
OS 



J.V. Alvarez and Felix Yndurain 

Departamento de Fisica de la Materia Condensada, 
Umversidad Autonoma de Madrid, 28049 Madrid, Spam 
(Dated: February 1, 2008) 

We discuss the fixed-point Hamiltonian and the spectrum of excitations of a quasi-bidimensional 
electronic system supporting simultaneously antiferromamagnetic ordering and superconductivity. 
The coexistence of these two order parameters in a single phase is possible because the magnetic 
order is linked to the formation of a metallic spin density wave, and its order parameter is not 
associated to a spectral gap but to an energy shift of the paramagnetic bands. This peculiarity 
entails several distinct features in the phase diagram and the spectral properties of the model, which 
may have been observed in CeRhlns. Apart from the coexistence, we find an abrupt suppression 
of the spin density wave when the superconducting and magnetic ordering temperatures are equal. 
The divergence of the cyclotron mass extracted from de Haas-van Alphen experiments is also 
analyzed in the same framework. 
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The interplay between magnetism and superconduc- 
tivity is a recurrent area of research in condensed matter 
physics. This interest has being activated in the last 
years due to the experimental findings of their coexis- 
tence in materials based on Ce, particularly in the 1-1- 
5 CeMIns family [H 0|- A prominent member of this 
family is CeRhlns, which grows in tetragonal form, al- 
ternating Celna and Rhln2 planes along the c crystallo- 
graphic axis. The structural anisotropy, induces quasi- 
bidimensionality in the electronic bonding and the Fermi 
surface, as evidenced in a series of de Haas- van Alphen 
measurements and band structure calculations 
At ambient pressure CeRhlns becomes an antiferromag- 
net (AFM) below Tn ~ 3.8K Q, with a small staggered 
magnetization aligned in the ab plane. Within the stan- 
dard Doniach's Kondo lattice paradigm, applying pres- 
sure in a weak AFM heavy-fermion system opens a route 
to very interesting effects. As the pressure increases this 
theoretical scenario predicts: i) A reduction of Tn due 
to Kondo compensation, ii) The eventual suppression of 
the AFM order in a quantum critical point (QCP), which 
alike other heavy-fermion compounds, would be respon- 
sible of the anomalies in the metallic phase, iii) The set- 
ting of unconventional superconductivity. However, all 
calorimetric NQR |7|], transport Q and susceptibil- 
ity 0] measurements provide a consistent picture for the 
pressure-temperature phase diagram (presented schemat- 
ically in Fig. [1]) in conflict with the aforementioned the- 
oretical scenario. Surprisingly, Tjv first increases with 
pressure and it only starts to decrease for pressures higher 
than 0.7GPa. Superconductivity shows up before T/v has 
gone down to zero, i.e. AFM and SC coexist. Finally, 
AFM disappears abruptly at Pc = 1.9GPa exactly when 
T/v = Tsci in a first order transition and before a QCP 
could have taken place. Despite the lack of a QCP in the 
pressure-temperature phase diagram, the metallic phase 
still might be understood in the framework of a quan- 
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FIG. 1: Schematic experimental pressure-temperature phase 
diagram of CeRhlns at zero-magnetic field adapted from ref- 
erences i, i and [ni. 



tum criticality if, changing another experimental knob, 
one could find a QCP nearby in the phase diagram. The 
natural choice is using a magnetic field to quench the 
superconductivity, and in that way, continue the point 
Pc{H = 0) into a line of first order transitions down to 
zero temperature, ending with a QCP at H*2. A new 
surprise appeared on this type of experiments 0, [lo| . 
For pressures higher than the AFM reenters applying 
an in-plane magnetic field Hrn< Hc2- Besides, trans- 
port measurements suggest [l|, Q that magnetic order in 
CeRhlns may not be associated to a gap in the single- 
particle spectrum. Actually, the resistivity as a function 
of the temperature does not show the conventional mini- 
mum characteristic of a metal-insulator transition at any 
temperature. The small anomaly observed in the resistiv- 
ity close to Tn seems related to a change on the scattering 
mechanism when the AFM sets in. 
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In vivid contrast with quantum critical and quasi-one- 
dimensional systems, the understanding of the individual 
AFM or SC states rely on simple but accurate mean-field 
theories. In this Letter, we propose that a basic compre- 
hension of the microscopic coexistence of AFM and SC 
and the first-order transition between these two phases 
can also be achieved within a mean-field scenario. We dis- 
cuss this phenomenon at the microscopic level in terms 
of a quasi-bidimensional model of interacting electrons 
proposed by one of us 12, [l^, We will enumerate 
some implications of that model, discussing to which ex- 
tent can be related to the phenomenology observed in 
CeRhlns. 

The model Hamiltonian contains the terms that nat- 
urally establish superconductivity and antiferromag- 
netism. 
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Where U is a Hubbard on-site repulsion, and Vkq the 
effective interaction in the Cooper channel, which we will 
take to be attractive. Only s-wave pairing will be con- 
sidered throughout this work. Superconductivity with 
an order parameter having a symmetry related to Vkq 
is favored by Hy but the presence of a Hubbard repul- 
sion establish a competition, which in terms of the gap 
is given by A = A2 — Ai, such that: 



A2 = Vj2{clJ-k-.) 

k 



(5) 
(6) 



where the prime restricts the summation to states with 
an energy (measured from the FL) smaller than a cut-off 
energy Ec. For those states we assume a very weak mo- 
mentum dependence of Vkq. In the absence of magnetic 
order and setting a constant density of states at the FL 
(i.e. far from a logarithmic divergence) , the competition 
between terms favoring and disfavoring the SC is clearly 
shown in the McMillan-like formula for the critical tem- 
perature. 
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FIG. 2: Calculated electronic structure of paramagnetic 
CeRhlns. (a) Density of states and (b) band structure. The 
arrow indicates a saddle point singularity at the X point of 
the Brillouin zone. The inset shows the variation of the sad- 
dle point energy along the z R-X-R direction (it crosses the 
Fermi Level at approximately the (vr/a, 0, O.Gtt/c) point). 



An essential ingredient of the model is the quasipar- 
ticle dispersion relation e(k) in the paramagnetic phase 
taken to set the Fermi level (FL) very close to the sad- 
dle point of a 2D system of itinerant electrons. Based 
on first principles calculations. Hall et al.Q find at the 
FL a sharp peak in the Density of States characteris- 
tic of a 2D van Hove logarithmic singularity. We have 
also performed a full first principles calculation of the 
electronic structure of tetragonal paramagnetic CeRhIn5 
using the VASP code 15|. The generalized gradient ap- 
proximation of Perdew et al.[l6| for the exchange and 
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correlation was adopted. The results for the experimen- 
tal lattice parameters are reported in Figure 2. In 
Figure 2(a) we observe, like Hall et al.Q, a sharp peak 
in the density of states near the Fermi level. In addi- 
tion, in Figure 2(b) we display the band structure at the 
vicinity of the Fermi level. We find a saddle point very 
close to the Fermi level whose dispersion along the z- 
direction (R-X-R in the standard notation) is about 0.3 
eV indicating the two-dimensional character of the sad- 
dle point singularity. Under these conditions we have, 
besides the Van Hove singularity in the density of states 
that favors superconductivity and magnetic order, an im- 
portant kinematic restriction, namely e(fc) = e{k + Q), 
where Q = (tt/o, tt/o) is the vector connecting two equiv- 



alent saddle points within the Brillouin zone. The above 
restriction determines the gapless nature of the SDW hke 
in the CDW case proposed by Rice and Scott [ist - 
The SDW order parameter is given by, 
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with = — 7-(T — 7- The resuhing quadratic Hamil- 
tonian is solved by a Bogoliubov transformation and the 
SC and SDW order parameters obtained selfconsistently. 
The Hamiltionian eigenvalues are obtained by solving: 
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The four solutions are Ei 

E.,E4, = E+ where E±{k) = ^e(fc)2 + A2 ± 7. No- 
tice that for a ID or 2D nested system e{k) = —e{k + Q) 
and therefore E±{k) = ±y/e{k)'^ + A^ +7^. 

The system is solved self-consistently for an effective 
half bandwidth at ambient pressure of Wq and the phase 
diagram for V = AWq, Ec = O.TWq , n = 0.92 electrons 
and U = 2.25Wo , U = 2.50Wo and U = 2.75Wo, is 
presented in Fig. [3l To simulate the effect of the pres- 
sure we have considered a linear variation of the band- 
width with the pressure and a bandwidth independent 
electron-electron interaction U. These assumptions have 
being found to be reasonable with a first principles calcu- 
lation using the SIESTA code [l9[ taking Ni as a bench- 
mark. Also, for simplicity, a constant V interaction and 
hole concentration independent of pressure are assumed. 
The results of our model agree qualitatively with the ex- 
perimental findings depicted in Fig. [T^nd Fig. 3(d). We 
indeed find a competition between SDW and SC but, as 
seen experimentally, they can coexist in a non negligi- 
ble region of the phase diagram. In addition, the SDW 
disappears abruptly when the two critical temperatures 
became equal, i.e., the SDW transition temperature can 
not be lower than the superconducting one. This numeri- 
cal result is similar to the analytical finding by Bilbro and 
McMillan (20| concerning superconductivity and marten- 
sitic transformation in A15 compounds. 

The proximity of the FL to a saddle point is an ingredi- 
ent of model ([T]) necessary for the formation of the metal- 
lic SDW. We ask ourselves whether in CeRhlns there is 
experimental evidence for such a proximity. A recent de 
Haas-van Alphen study j3] shows a divergence in the cy- 
clotron mass at a pressure Pc2 ~ 2.4GPa, accompanied 
with a change in the quasi-2D Fermi surface. Those ex- 
periments were performed for values of magnetic fields 
and pressures in which the system is AFM. Following 
Park et al. Pc2 is very close to the pressure at which T^r 
would extrapolate to zero in absence of SC (see upper 
inset in Fig. 2] ). The cyclotron mass in a 2D elec- 
tronic system nic = /27r{dA{E) /dE) where A is the 



_ 


' 1 ' 1 


(a)- 




- 






- 


_ • 
• 
• 

- • 
• 


,1,1 


, 1 


1 2 


3 4 


5 



(c)_ 



-' ll' I ',1 ' I ' I ' I ' I ' ' I ' I ' I '- 






1 ' 1 

(b): 

- 







2 3 


4 5 


' 1 ' 

-T(K) 


' 1 ' 1 

\, Exp.J 


mm m • # , 



WAV„ 



5 10 15 20 25 
P 



FIG. 3: Phase diagram (temperature versus pressure) ob- 
tained using the model Hamiltonian and a two dimensional 
band structure with the FL close to a saddle point. The vari- 
ables are in units of the half bandwidth Wo- Tn and Tsc 
stand for the SDW and SC critical temperatures respectively. 
Panels (a), (b) and (c) are the calculated results for the pa- 
rameters given in the text and U = 2.25Wo , U = 2.50VKo and 
U — 2.75W0 respectively. The inset in panel (c) indicates the 
metaUic DOS in the three (SDW, SDW-hSC and SC) different 
regimes. The shaded area indicates occupied levels. Panel (d) 
represents the experimental results of ref. 

area enclosed in the isoenergetic contour line E{k) — E. 
Close to the 2D Van Hove singularity one expects the 
cyclotron mass to diverge logarithmically. Actually, the 
precise functional form has been computed in Ref. 
and found to be 
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In a tetragonal crystal structure C and D are numbers 
nearly independent of the pressure and rrico is the cy- 
clotron mass at the bottom of the band. The divergence 
is driven by the denominator in the argument of the log- 
arithm which in our model is: Ep{'-f) — Eyhs ^ l{P) ^ 
Tn{P)- In other words, in model ([1]), rric is enhanced 
as the AFM vanishes because the FL in the SDW ap- 
proaches the saddle point at X (see the inset in Fig. [4]). 
To elucidate if the experimental results are compatible 
with this argument and with expression ^ , we have ex- 
tracted the pressure dependence of T^r, fitting the ex- 
perimental data in Ref (Q) (see upper inset in Fig. [4] 
to a cubic polynomial law in {Pc2 — P) in the range of 
pressures between P = 0.65 and P = 2.4 GPa. We do 
not attempt to justify physically this fit here, since our 
only goal is extracting an analytic expression for T]y{P) 
for the range of pressures of interest, to insert it in ([9]). 
The results are presented in the main panel of Fig. |4l 
Our model reproduces reasonably well these experimen- 
tal findings. 

Within this model, we expect an anomaly in the spe- 
cific heat at Tn, which is not associated to a SDW spec- 
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FIG. 4: Cyclotron mass {rric) as a function of pressure. Tri- 
angles are the experimental values from Ref. [4j] showing the 
nic divergence close to P ~ 2.4 GPa, where the magnetic 
order disappears. The solid line is a fit to the theoretical 
model. Close to a Van Hove singularity the cyclotron mass 
diverges logarithmically as the difference between Ep and the 
Evhs vanishes (see lower inset). According to the model this 
energy difference is proportional to Tat and its pressure de- 
pendence can be extracted, for instance, from Ref. 0](see 
upper inset). 

tral gap but to the entropy released when the magnetic 
order disappears. The electronic part of the specific heat 
is: 
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where i=3,4. The second term inside the parenthesis 
gives the SC anomaly at Tc and the third term gives 
the antiferromagnetic anomaly at T/y. The SC anomaly 
is much weaker in the coexisting phase, because the FL 
hes in a depression of the density of states (see central 
graph on the panel (c) of Fig. [3]) created by the under- 
lying SDW, while in the purely SC phase the FL is very 
close to a divergence in the density of states (right graph) . 
Remarkably, this enhancement has been also observed in 
calorimetric measurements on CeRhlus Q 

To summarize: beyond the detailed boundary shape 
in the phase diagram of CeRhlns we have identified two 
unequivocal, clear-cut features of the phenomenology of 
CeRhIn5 in the discussed model, namely; coexistence and 
abrupt disappearance of AFM when Tn = Tsc- The es- 
sential ingredient in the model is the metallic SDW, fa- 
vored for the proximity of Fermi level to a Van Hove 
logarithmic singularity in the density of states. The 
gapless nature of the SDW implies the lack of a metal- 
insulator transition at, or close to, Tj\j as shown by resis- 
tivity measurements. The kinematic conditions needed 



for the metallic SDW to be formed seem to be present in 
CeRhlns as shown by the logarithmic divergence of the 
cyclotron mass. 
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